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We investigate the electronic transport through two parahel double quantum dots coupled both 
capacitively and via a perpendicularly aligned charge qubit. The presence of the qubit leads to a 
modification of the coherent tunnel amplitudes of each double quantum dot. We study the influence 
of the qubit on the electronic steady-state currents through the system, the entanglement between 
the transport double quantum dots, and the back action on the charge qubit. We use a Born- 
Markov secular quantum master equation for the system. The obtained currents show signatures of 
the qubit. The stationary qubit state may be tuned and even rendered pure by applying suitable 
voltages. In the Coulomb diamonds it is also possible to stabilize pure entangled states of the 
transport double quantum dots. 



I. INTRODUCTION 

The capability to harness the potential power of 
quantum properties, such as superposition of states 
and entanglement, for information and communication 
technologies, has raised interest in possible experimental 
implementations over the recent years. Scalability, the 
feasibility of coherent control, and non-destructive read 
out of quantum states as well as robustness against 
decoherence are among the key features for such can- 
didates [ ]. Encouraged by the latest progress in their 
fabrication and manipulation, semiconductor quantum 
dots (QDs) have been proposed as possible candidates 
[: ] and have been intensively studied since then. 
Of special interest are double quantum dots (DQDs) [ ], 
which are used to model two- level qubit states [ - ]. It 
has been shown that coherent control [ - ] and read 
out [x^-^.^] are achievable in these systems. Together, 
this allows us to implement electronically accessible 
quantum gates based on quantum dots [14-16]. Even 
entanglement, which is crucial for quantum computation, 
can be produced [ ], manipulated [18] and detected 
[19, 20] in quantum-dot setups. In consequence, quan- 
tum dots provide a promising candidate for a quantum 
computation architecture. 

Nowadays, quantum dots are produced on a large 
scale, for example, via lithographic methods [21, 22], 
self-assembled growth [ ], or by depletion of two- 
dimensional (2D) electron gases in semiconductors [ ], 
to name but a few. Unfortunately, all these fabrica- 
tion methods can not completely exclude unwanted 
impurities that might destroy the desired properties 
of double quantum dots. The effects of impurities on 
coupled quantum-dot systems have been recently studied 
theoretically using molecular orbital and configuration 
interaction methods [ ]. Furthermore, the possibility 
of screening charge impurities by using multi-electron 
QDs were theoretically investigated [ ]. In addition, 
spin impurities have been observed experimentally via 
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transport spectroscopy in a carbon nanotube DQD [^'^]. 
The main goal of this paper is to analyze the effect of 
impurities on the transport and system characteristics 
of otherwise ideal DQDs. In contrast to the previously 
mentioned methods, we are using an effective rate 
equation approach derived from a Lindblad master 
equation. This ansatz yields a practical means to theo- 
retically derive the transport properties of such systems, 
in particular, the transport spectra. To this end, we 
investigate a setup of two parallel DQDs, which give rise 
to entangled two-particle states, and a central impurity 
that disturbs the intrinsic DQD tunneling amplitudes. 
Already without impurities, these systems may exhibit 
interesting non-standard fluctuation relations [28]. 
In order to keep the calculations simple, we model 
the impurity as a charge qubit (CQB), which perturbs 
the DQD system depending of its internal state. The 
respective system Hamiltonian is introduced in Sec. II A. 
Subsequently, using the master-equation approach 
described in Sec. II B, we analyze its influence on the 
electronic transport spectrum through the DQDs both 
numerically and theoretically within Sec. III. Finally, we 
study the purity of the system in Sec. IV and investigate 
the effect of the CQB on the entanglement of the 
two-electron states in Sec. V. 



II. MODEL 

We consider an experimental setup as depicted in 
Fig. 1, which consists of two parallel DQDs that are cou- 
pled both capacitively and via a perpendicularly aligned 
CQB impurity. Furthermore, these DQDs are coupled to 
leads that act as reservoirs for electrons and, thus, al- 
low for electronic transport when their parameters, i. e. , 
temperature and/or chemical potential, are chosen differ- 
ently. To simplify the calculations, we explicitly assume 
spin-polarized electronic leads and neglect the spin de- 
gree of freedom throughout this paper. This implies that 
spin-selection effects as observed in Ref. [ ] do not mat- 
ter. 
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Figure 1. (Color online) Two double quantum dots (DQDs) 
1 and 2 with electronic on-site energies e, coherent tunnel- 
ing amplitude (solid lines) Tc, and on-site Coulomb interac- 
tion (dashed lines) U\\ interact capacitively via perpendicular 
and diagonal Coulomb interactions U± and Ux, respectively. 
Transport is enabled by tunnel coupling the quantum dots 
(i, a) to adjacent fermionic reservoirs of chemical potential 
fii^a with the tunneling rate Ti^a- A charge qubit (CQB) 
with detuning Q and tunneling amplitude J modifies the tun- 
nel amplitudes of both DQDs by G (wavy lines) depending 
on the position of the CQB electron. 



Hamiltonian 



The Hamiltonian of the full system can be decomposed 
into H = Hs + Hb + ^sb with a system Hamiltonian 
a bath Hamiltonian i^B and a Hamiltonian Hsb de- 
scribing the coupling between system and bath. In the 
considered setup, the electronic transport through each 
DQD is altered by the state of a nearby impurity which 
is modeled as a CQB. In particular, the current through 
the DQD closest to the charge of the CQB is suppressed 
due to Coulomb repulsion. 

We define the electronic occupation number operator 
fia^ = d] ^di^a where the operators di^a and d] ^ annihi- 
late and create electrons in the quantum dot a G {R^ L} 
in channel i G {1,2}. For convenience, we use here and 
in the following discussion the notation z = 1 for the up- 
per and z = 2 for the lower transport channel in Fig. 1, 
whereas the labels L, R denote the left and right quantum 
dots and electronic baths, respectively. Furthermore, we 
introduce the operators q, c] which annihilate and cre- 
ate electrons in the CQB. Since the CQB is equivalent 
to a two-level system, we can express the corresponding 
operators in terms of Pauli matrices a 
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The first two lines in Eq. (1) describe the isolated DQD 
properties. Specifically, the first term accounts for the 
homogeneous on-site energies e of the DQDs. The sec- 
ond term describes the intrinsic tunneling processes in 
the transport DQDs, which are proportional to the co- 
herent tunneling amplitude Tc > 0. The term in the 
second line stems from the Coulomb interaction within 
each DQD that is parametrized by the Coulomb interac- 
tion strength U\\. In the third line of Eq. (1), we include 
the Coulomb interactions in-between the parallel DQDs 
that arise from electrons on the same side with strength 
U±_ and from electrons on different sides with strength 
Ux. Finally, the last two lines of Eq. (1) describe the 
CQB with detuning Vt and coherent tunneling amplitude 
J, and its influence on the transport DQDs: The intrin- 
sic processes of the DQDs are modified by the parameter 
G in the presence of an electron in the respective dot of 
the CQB. Due to Coulomb repulsion, it is reasonable to 
assume that the presence of an electron in the CQB leads 
to a suppression of the tunneling amplitude Tc. Thus, we 
restrict the parameter G to negative values G G [— Tc, 0]. 
The electronic baths are modeled as free electron gases 
of spin-polarized electrons. Hence, the Hamiltonian for 
bath (z,a) is given by 



Hr 



^i,a,k Qi a k Q'iyOi-,k-) 



(2) 



with operators qj ^ ^ and qi^a.k creating and annihilating 
an electron with momentum k and energy Vi^a.k in lead 
(z, a). The system-bath interaction Hamiltonian reads as 



^SB 



so 



H 



(3) 



and cFx = cjc2 + c^ci. Using these abbreviations yields a 



where the tunneling amplitude of an electron hopping 
from the lead (i,a) into the system or vice versa is pro- 
portional to t*^ ^ and ti^a,k^ respectively. 



B. Liouvillian 

Assuming that the leads are in thermal equilibrium 
and the coupling between the leads and the system is 
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weak, we can use the Born-Markov secular approxima- 
tion (BMS) [ ]. Starting from the von Neumann equa- 
tion p — i2 [p] , this formahsm ahows us to extract a 
quantum master equation that assumes the form of a 
rate equation for the reduced system density matrix p in 
the system energy eigenbasis for non-degenerate energy 
eigenvalues. Within this ansatz, the Liouville super op- 
erator C is parametrized by the Markovian system-bath 
tunneling rates Vi^a.{y^) = 27r^^\ti^a,k\'^ - ^i,a,k) 
and the Fermi functions of the leads fi^ai^k) = 

^^B{qj^a,k^^,cx,k Pb} = 1/ {exp [/3 - Pa)] + !}• 

Here, we introduce the inverse temperature /3 = 
1/ (/cbT), which we assume identical for each lead, and 
the chemical potential pa foi" the left or right leads. For 
sequential electronic tunneling, we can uniquely identify 
the jump terms in the master equation which enables one 
to convert it into a conditional master equation. Consid- 
ering the number n of electrons tunneled via one lead 
(z, a) into or out of a system with non-degenerate energy 
eigenvalues, the conditional master equation reads as 

p^^) = Co p(^) + + C- . (4) 

Here, the super-operators £o, >C^, and jC~ are acting 
on the reduced system density matrix. The operators 
/ jC~ describe electronic jumps into/from the moni- 
tored reservoir from/into the system, respectively. The 
operator Cq describes the internal dynamics of the sys- 
tem and electronic jumps between the system and the 
remaining baths. This n-resolved master equation can 
also be established using virtual detectors as bookkeep- 
ing operators [ ]. 

Subsequently, we perform a Fourier transformation 
PiXa^t) = exp(i nxa), which introduces a 

counting field Xa for the respective lead a. Hence, the 
Liouville super operator for lead a becomes a function of 
this counting field: 

C{xa) = Co^C-^e^'^-^C-e-'^-. (5) 

This procedure can be applied for each of the four elec- 
tronic baths leading to a corresponding number of count- 
ing fields. However, due to charge conservation and the 
neglect of short-time dynamics, we just need to consider 
one counting field for each transport channel. Hence, 
without loss of generality, we consider the leads (l^L) 
and (2, L) and denote the corresponding counting fields 
as xi and X2- In the following calculations, we focus on 
the steady-state currents through the system. For lead 
a, this current [30] is defined by the relation 

= -iTr{ax„£(xi,X2)|^,=^,=o 4' ^^'^ 

where the steady-state reduced density matrix p is deter- 
mined by = £(0, 0) p. Finally, we note that throughout 
this paper, we utilize the fiat-band limit with energy- 
independent tunneling rates Ti^^ {^) = ^i,a- 



III. TRANSPORT CHARACTERISTICS 

Within the following sections, we will apply the master 
equation formalism derived above to the Hamiltonian (1) 
and analyze its steady-state properties both analytically 
and numerically. 

A. Full Transport Characteristics 

First, we analyze the steady-state transport spectrum 
of the full Hamiltonian (1). To this end, we evaluate 
Eq. (6) numerically. This requires to set up the full 
32 X 32 Liouvillian (see Appendix A). The transport spec- 
trum is obtained by calculating the steady-state current 
for varying gate voltage Vcate = £ which shifts the en- 
ergy levels of all transport DQDs according to Eq. (1), 
and varying external bias voltage ^Bias- The bias volt- 
age enters the Liouville super operator through the Fermi 
functions via the chemical potentials. For convenience, 
we assume symmetric chemical potentials /iL = '^Bias/2 
and /iR = — FBias/2 for both transport channels in all 
further calculation. As an example, we show the differ- 
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Figure 2. (Color online) Unitless differential steady-state cur- 
rents through channel 1 (left) and channel 2 (right). Most 
prominent are the CD structures with a different number of 
electrons (numbers) loaded into the system. The dashed cir- 
cles indicate the presence of gaps between the CDs and the 
dashed vertical lines correspond to Fig. 3. The other param- 
eters are fixed to: /3 = 100/Tc, U\\ = 20Tc, U± = Ux = 10 Tc, 
J = Q = Ti^c. = Tc and G = -Tc. 

ential transport spectra for both transport channels in 
Fig. 2. The largest differential currents dl /dVBias ^ire 
observed at the edges of the Coulomb diamonds (CDs). 
Within these regions, the steady-state current vanishes 
due to Coulomb blockade effects that prevent electronic 
transport through the system [ ]. For a suitable choice 
of gate and bias voltages, this blockade can be overcome 
leading to finite steady-state currents. The Coulomb di- 
amonds can be associated with specific electronic sub- 
spaces reaching from the vacuum state with zero elec- 
trons up to the maximal number of four electrons, which 
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are indicated in Fig. 2. Note that distinct CDs can only 
be observed for low temperatures. 



B. High-Bias Currents 

Unfortunately, the Hamiltonian (1) is complicated to 
analyze analytically. Thus, we first consider the case 
of a high bias voltage Vbias ^ 0, where only the sys- 
tem transition energies from the zero-to-one electron sub- 
space lie within the transport window. In this ultra- 
strong Coulomb blockade (USCB) regime, we make use 
of Eq. (6) and find that the high-bias steady-state current 
through the transport channel i becomes 



^USCB — 



2 (r.,Lr^ 



(7) 



i,R 



where i labels the opposite transport channel. 
For comparison, we additionally calculate the high-bias 
steady-state current for the strong Coulomb blockade 
(SCB) regime where in each DQD at most one electron is 
allowed. In this high-bias limit, where all transition en- 
ergies from the zero-to-one and one-to-two electron sub- 
space lie within the transport window, the steady-state 
current through transport channel i is of the form 



^SCB — 



2F, 



i,R 



(8) 



This result corresponds to the high-bias steady-state 
current one obtains for sequential electronic transport 
through a two-level system [3i] in the SCB regime. From 
Eq. (7), we see that in the USCB regime the current 
^USCB explicitly depends on the tunneling rates of the 
opposite channel i. This is to be expected for this config- 
uration with dynamical channel blockade [3"^, 88]. This 
intermediate coupling is lifted, for example, if the DQD 
i is almost immediately reloaded from the left lead, i. e. , 
^i,L ^ ^i,R- Then, the system tends to be always occu- 
pied by an electron: The steady-state current in channel 
i becomes proportional to the coupling to the right lead 
Ti^R^ whereas the current through the other transport 
channel vanishes. In the opposite case, when the DQD i 
is not refilled from the left lead, i. e. , Ti L 0, while the 
other tunneling rates remain non- vanishing, the steady- 
state current through channel i vanishes and the other 
current takes on the form of Eq. (8). Considering a sim- 
ilar configuration where the electrons almost instantly 
leave the transport channel i via the right lead, i. e. , 
^i,R ^ ^i,L^ we find that the channel coupling is only 
partially lifted. In fact, in this limit, channel i decou- 
ples and the respective steady-state current takes on the 
form of Eq. (8). However, the steady-state current i is 
still proportional to tunnel couplings of both transport 
channels. Finally, if the couplings to the left leads and 
the couplings to the right leads are the same for both 
channels, i. e. , F^^q, Fq,, the steady-state currents are 
the same for each channel. 
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Figure 3. (Color online) Plot of the steady-state currents 
for both transport channels 1 (dashed) and 2 (solid). The 
plot corresponds to the dashed vertical line in Fig. 2 at 
Voate = 0.8 Tc. We indicated the high-bias currents in the 
ultra-strong Coulomb blockade (/uses) and strong Coulomb 
blockade (/scb) regimes. The insets compare the steady- 
state current of channel 2 (solid) with the results for par- 
allel DQDs without impurity (dashed) and a single DQD 
coupled to a CQB (dotted). The system parameters are 
fixed to: /3 = 100/Tc, U\\ = 20 Tc, U± = Ux = 10 Tc, 
J = Q = Ti^c. = Tc and G = -Tc 



This demonstrates that in the high-bias regime, the 
steady-state currents are not sensitive to the asymme- 
try induced by the CQB. Note that due to electron- hole 
symmetry, we find analogous results to Eq. (7) if the sys- 
tem transition energies from the three-to-four electron 
subspace lie within the transport window. 



C. Current Anti-Correlation 

In this section, we further analyze the transport prop- 
erties of the system in the ultra-strong Coulomb block- 
ade regime for small bias voltages. Thus, we focus on 
the region in the vicinity of the edge of the 0-electron 
CD of the transport spectra in Fig. 2. A first interest- 
ing feature in this region is the occurrence of negative 
differential conductance, which indicates blocking effects 
in both transport channels. We demonstrate in the inset 
in Fig. 3 that this is an intrinsic feature resulting from 
the CQB impurity as negative differential conductance 
is not present in parallel DQDs without impurity. Ef- 
fectively, the asymmetry induced by the impurity leads 
to the fact that the two possible transport channels be- 
come accessible at different gate voltages. Furthermore, 
since the total number of electrons in the system is con- 
strained by Coulomb interactions, as soon as the sec- 
ond transport channel becomes available, the current in 
the first channel decreases, leading to a negative differ- 
ential conductance. Consistently, this property is most 



prominent in the Coulomb blockade regime. Therefore, 
by comparing the differential currents for the two trans- 
port channels, we find that this feature is accompanied 
by an anti-correlation between the steady-state currents 
in the different transport channels. 



D. Coulomb Diamond gaps 

The most striking difference between the plots in Fig. 2 
is the fact that some Coulomb diamonds do not close, 
as indicated by the dashed circles. We expect this ef- 
fect to occur for coupled DQDs with asymmetries in the 
quantum-dot energies or in their tunneling amplitudes. 
Because we explicitly exclude this kind of asymmetry in 
the transport DQDs, this effect clearly suggests a block- 
ing induced by a finite energy barrier at zero bias voltage 
that stems from the presence of the CQB. In which chan- 
nel the gap appears depends on the sign of the detuning. 
Since all changes in the steady-state currents are associ- 
ated with resonances of the transition energies between 
eigenstates of the system Hamiltonian, it is possible to 
calculate the position of the lines shown in Fig. 3 if the 
eigenvalues of the respective Hamiltonian are known. 
Hence, a diagonalization of Eq. (1) allows one to calcu- 
late the position and width of the gap. For example, the 
gate voltage with minimum gap between the 0-electron 
and 1-electron CD is given by 

^cTe = ^ (2 G + 4 Te + x/4J2 + (G-!^)2 

-2 J^ + 9? + a/4 J2 + (G + f2)2) . (9) 

Subsequently, we derive the other parameters of the gap 
such as the upper bias voltage at this point, which reads 



Trlp'^) 



T^upper 
'^Bias 



'J2 



G + Vl 



■ (10) 



Due to the choice of our parameters, the gap is symmetric 
with respect to T^Bias = and thus the lower bias voltage 
satisfies V^;^^' 
the gap 



"^^s^^' which yields for the width of 



^Gap 



V4 J2 + (G - - v/4 J2 + (G + ^f. (11) 



In an analogous way, we can determine the point of con- 
tact of these Coulomb diamonds. We find that this posi- 
tion is given by the relation 



^Gate 



G 

~2 




From Eq. (11), we see that for the considered experimen- 
tal setup, the appearance of the gap results from both, 
the detuning Q of the CQB and the modification G of the 
intrinsic tunnel amplitudes. If one of these quantities is 
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Figure 4. (Color online) Density plot of the purity of the 
system density matrix for low temperatures /3 = 100/Tc. The 
picture also shows the respective normalized eigenvectors of 
indicated in the CDs by So ,Si, &, S3, and S4, which 
correspond to Eq. (13). The system parameters are fixed to 
U\\ = 20Tc, U± = Ux = lOTc, J = Q = Fi^c. = Tc and 



zero, the gap vanishes. On the contrary, the width of 
the gap is completely independent of the intrinsic tunnel 
amplitudes Tc of the transport DQDs. Moreover, we find 
that the gap not only vanishes as G or 1] approaches zero, 
but also if J becomes very large. This behavior can be 
understood since the rapidly oscillating CQB on average 
affects both transport channels in the same way and the 
energy barrier vanishes. Investigating the position of the 
gap according to Eq. (9) yields a more complicated be- 
havior in dependence of the CQB parameter J. As the 
hopping amplitude J is increased from zero, the gap is 
shifted from an initial finite value to higher gate voltages 
up to a maximal value. A further increase of J results in 
a shift of the minimum gap to lower gate voltages, which 
in the limit J ^ 00 becomes a constant that equals the 
position for J = 0. Therefore, we conclude that the sys- 
tem is most sensitive to the CQB for a small tunneling 
amplitude J. 

In general, we note that measuring in an experiment the 
quantities described by Eqs. (9)-(12) allows us to calcu- 
late the tunneling amplitude J and the detuning Q of the 
CQB as well as the modification G of the tunnel ampli- 
tudes Tc of the DQDs. 



IV. PREPARATION OF PURE STATES 

The SCB regime is particularly interesting because it 
includes configurations where the transport DQDs can 
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also be treated as CQBs, giving rise to an effective system 
composed of three coupled qubits. These coupled qubits 
represent possible candidates for basic constituents of 
quantum information technologies that use entangled 
qubit states for calculations or communication. In order 
to better understand the eigenstates of the considered 
system in the long-time limit, we analyze the purity 
of its reduced steady-state system density matrix and 
investigate some interesting expectation values. 



A. Purity of the Full System 

In Fig. 4, we plot the purity Tr {p^} of the reduced sys- 
tem density matrix. We find that in the interior of the 
Coulomb diamonds and for low temperatures, the sys- 
tem enters a pure energy eigenstate with the lowest pos- 
sible energy. This behavior causes the vanishing steady- 
state currents in this region. In general, these eigenstates, 
which are obtained from a diagonalization (see Appendix 
A) of the Hamiltonian (1), are entangled states. However, 
for the special choice of Coulomb interactions Uj_ = f/x, 
we find that the eigenstates become separable. There- 
fore, the eigenstates Si the system enters in the CDs of 
Fig. 4 have the simple representation 

5o=|O,O)0(ao|;)-6o|t)), 

^i=(|L,O)-|i?,O))0(ai ID), 

S2 ={\L)- \R)) - \R)) (a2 |t) - , 

^3 = {\L, LR) - \R, LR)) (as \i) - bs |t)) , 

^4 = \LR, LR) (a4 ID - 64 It)) , (13) 

where the structure of the kets is defined in the local 
basis as |DQD]^, DQD2, CQB) and {ai^bi} G R are the 
respective normalized coefficients. 

If the temperature is raised, the region within a CD 
where purity is reached shrinks in favor of a mixture 
of states with the same number of electrons. Outside 
of the Coulomb diamonds, eigenstates belonging to dif- 
ferent electron subspaces mix together, which allows for 
electronic transport and a finite current. 
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Figure 5. (Color online) State analysis of the charge qubit 
for low temperatures /3 = 100/Tc. The left picture shows the 
modulus of the Bloch vector and the right picture shows its az 
component. The system parameters are fixed to U\\ = 20 Tc, 
U^ = Ux = lOTc, P. a = J = Q = Tc and G = -Tc. 



one. Because in our model there is no ay component for 
the CQB, the modulus of the Bloch vector is defined by 

II {^) II = 



(^x)^ + (<^z)^- We find that it is possible to 
tune the CQB to a pure qubit state for suitable combi- 
nations of gate and bias voltage (see Fig. 5). The regions 
with large CQB purity coincide with the Coulomb dia- 
monds indicating that the CQB is not entangled with the 
rest of the system. Furthermore, there are well defined 
regions outside of the Coulomb diamonds where the CQB 
behaves completely classical, indicated by the vanishing 
of the Bloch vector modulus. However, as depicted in 
Fig. 5 the eigenstates that the CQB can take on within 
the Coulomb diamonds, do not correspond to localized 
electron states as the a^-component never takes on the 
value ±1 in these regions. These properties remain qual- 
itatively unchanged if other parameters are considered. 
In general, we find that the purity properties are sensi- 
tive to thermal fiuctuations and thus the stability of the 
corresponding regions will decrease as the temperature 
increases. This effect is strongest for the 2-electron CD 
and weakest for the 4-electron and 0-electron CDs. 



B. Purity of the CQB 

The observation that there exist regions where the 
whole system is in a pure state leads to the question if the 
same is true for the CQB. A pure system density matrix 
together with a mixed CQB density matrix would imply 
entanglement between the transport DQDs and the CQB. 
In order to answer this question, we calculate the station- 
ary expectation values of the CQB occupation operators 
according to ((jj) = Trja^ p}. Since the CQB represents 
a two-level system, these expectation values can be in- 
terpreted as the components of a Bloch vector. A pure 
CQB state corresponds to a Bloch vector with modulus 



V. ENTANGLEMENT 

In contrast to the separable eigenstates in Eq. (13), we 
can also change the Hamiltonian parameters to stabilize 
entangled states in the 2-electron CD. Thus, motivated 
by investigations of the entanglement of the system (1) 
without the CQB impurity [18, 34], we explore the effect 
of the presence of the CQB on the entanglement of the 
transport channels within the following section. 
In order to qualitatively and quantitatively determine 
the entanglement between the two transport DQDs, we 
project the steady-state density matrix onto the two- 
qubit subspace where exactly one electron is present 
in each transport DQD. Subsequently, tracing out the 
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Figure 6. (Color online) The left picture shows the concur- 
rence of the four energy eigenstates in the qubit sector of 
Hamiltonian (1) as a function of the ratio a of the Coulomb 
interaction strengths. For comparison the inset shows the 
concurrence for the system without CQB impurity. The right 
plot shows the generalized transport concurrence with a = 0.7 
for the full transport spectrum. All parameters are fixed to 
13 = 100/Tc, U\\ = 20 Tc, U± = 10 Tc, r^,c, = J = 11 = Tc and 
G = -Tc. 



CQB degrees of freedom yields an effective 4x4 matrix 
p2 = Trrest {p} for the two coupled qubits represented 
by the transport DQDs. Here, the trace over "rest" in- 
cludes the 0-, 1-, 3-, and 4-electron subspaces as well as 
the 2-electron states corresponding to a doubly occupied 
transport DQD. For this effective system of two coupled 
qubits, there exists a well-known entanglement measure 

[35]: The concurrence C = max 0, \/A7 — X]t=2 V% 

where are the eigenvalues of P2{o'y ^ ^y)p2{^y ^ ^y) 
arranged in decreasing order, i. e. , A^+i < A^. 



A. Eigenstate Concurrence 

For the system of two coupled transport DQDs it is 
known that the strength of entanglement depends on the 
strength of the on-site Coulomb interaction [ ]. Addi- 
tionally, the entanglement of two parallel DQDs is also 
induced by asymmetric interactions. Hence, this prop- 
erty strongly depends on the choice of Coulomb interac- 
tions between the quantum dots. In our numerical eval- 
uation in Fig. 2 we explicitly assume that the Coulomb 
interaction between parallel dots has the same strength 
as the diagonal interactions, i. e. , /7x = /7^. However, 
with this choice, entanglement can hardly be achieved. 
In order to investigate the dependence of the entangle- 
ment of the pure states on the interaction strength, we 
introduce the ratio a = Ux/U^ between parallel and di- 
agonal Coulomb interaction in the system Hamiltonian 
(1). First, investigating the concurrence of the pure qubit 
eigenstates for parallel DQDs with and without CQB in 
dependence of this ratio a yields the results presented 
the left picture in Fig. 6. In general we find that for the 
highly asymmetric case of a = the concurrence for par- 
allel DQDs with CQB impurity is maximal and reaches 



almost 1. With increasing ratio a the concurrence de- 
creases until it completely vanishes at a = 1, the point 
of maximum symmetry. In contrast, the two Bell states 
1/V2 {\R,R) - \L,L)) and - \L,R)) corre- 

sponding to the straight line in the inset in Fig. 6 have 
a constant concurrence of 1 over almost the full range 
of a. Hence, these eigenstates are maximally entangled 
states except near the point a = 1 where the concur- 
rence vanishes discontinuously. Note that always two of 
the four qubit states show exactly the same concurrence. 
Thus, both the solid and dashed lines in the inset cor- 
respond to two eigenstates each. Comparing with the 
concurrence for qubit eigenstates of Eq. (1), we find that 
the Bell states are destroyed due to the influence of the 
tunneling amplitude J and the detuning Vt of the CQB. 
This result is in agreement with recent theoretical work 
using a configuration interaction method to analyze im- 
purity effects on quantum bits [ ]. In addition, we ob- 
serve that in general the concurrence of non-Bell states 
is raised for all values of a. This effect mainly results 
from the modification parameter G. The concurrence is 
maximally enhanced for G = — Tc and approaches the 
value of the unperturbed system as G approaches zero. 



B. Generalized Transport Concurrence 

Next, we investigate the concurrence for the whole 
transport spectrum, which includes mixed states as well. 
Since both cases with a = and 1 are experimentally 
hard to achieve and also represent rather special config- 
urations, we choose a = 0.7 for numerical investigations. 
From the definition of the concurrence, it is obvious that 
this quantity is linear in the norm of the reduced density 
matrix p2. Hence, if the reduced density matrix is not 
renormalized after tracing out the CQB, the concurrence 
is only exact if evaluated in the transport qubit sector 
where a single electron is on each transport DQD. 
Otherwise, the calculation yields a product of the exact 
concurrence multiplied by the probability for being in a 
transport qubit state. Due to this special property we 
will in the following use the non-renormalized concur- 
rence G to characterize the entanglement of the DQDs 
in the full gate and bias regimes. In the right picture of 
Fig. 6, we plot the resulting concurrence versus gate and 
bias voltage. 

We observe that the presence of the CQB enlarges the 
concurrence within the 2-electron Coulomb diamond 
and slightly decreases the concurrence of the mixed 
states outside of this Coulomb diamond. Moreover, 
we see that in the exterior of the 2-electron Coulomb 
diamond, the concurrence vanishes almost everywhere 
except for the regions associated with transitions from 
the 1- or 3-electron to the 2-electron subspace. Here, 
the transport qubit eigenstates enter the steady-state 
reduced density matrix contributing their high con- 
currence to the mixture of states. However, since we 
do not re-normalize the effective two-qubit density 
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matrix, this high concurrence gets multiphed by the 
2-electron fraction of the mixed state which corresponds 
to the probability to be in such a 2-particle state. 
Going deeper into the 2-electron subspace allows for 
more 2-electron eigenstates to mix in the steady state 
reduced density matrix and hence rapidly reduces 
the respective entanglement. In an experiment the 
DQD, entanglement could be measured for example via 
construction of a Bell inequality as suggested in Ref. [18]. 



VI. SUMMARY 

We studied the electronic transport properties of a sys- 
tem of two parallel DQDs which are subject to a CQB im- 
purity that affects the transport through these DQDs. In 
particular, we investigated its effect on the steady-state 
currents in the strong and ultra-strong Coulomb block- 
ade regimes. We find that the CQB detuning gives rise to 
an asymmetry that prevents the closing of Coulomb dia- 
monds. Width as well as the position of the gap depend 
on the CQB parameters. In the ultra-strong Coulomb 
blockade regime we were able to extract analytic results 
for both, the steady-state current in the high-bias regime 
as well as the width and position of the gap in depen- 
dence of the CQB parameters. Moreover, we analyzed 
the purity of the reduced system density matrix and 
found that the purity of the system is preserved. In addi- 
tion, the back-action on the CQB allows us to even render 
its eigenstates pure. Finally, we analyzed the impact of 
the CQB on the entanglement between the two transport 
channels. Here, we observe on the one hand a destruc- 
tion of Bell states due to the CQB tunneling amplitude 
J and detuning Vt. On the other hand, we find an in- 
crease of entanglement of the remaining eigenstates due 
to the modification parameter G. In the exterior of the 
Coulomb diamond, entanglement is slightly decreased by 
the CQB. 

Financial support by the DFG (SFB 910, SCHA 1646/2- 
1) is gratefully acknowledged. We would also like to 
thank S. Kohler for helpful discussions. 



Appendix A: Explicit Liouvillian 

The system Hamiltonian defined in Eq. (1) has a block 
structure in the local basis {|DQD]^, DQD2, CQB)} where 
each block is associated with a fixed number of electrons 
in the transport DQDs. In order to derive the explicit Li- 
ouvillian we first diagonalize the full system Hamiltonian 

(1) 



Hs |n,a) = £„,a \n,a) . 



(Al) 



Here n G {0,1,2,3,4} labels the number of electrons 
in the transport channels. The second index a G 
{1, . . . , a^^^} with a^^"" = a^^^ = 2, = a^^^ = 8 

and af^^ = 12 labels the different eigenstates belonging 



to the same electronic subspace n. For convenience we 
first sort the basis vectors {|n, a)} of the eigenbasis of 
^^s by their increasing electron number n. Subsequently, 
the eigenstates belonging to the same number n of elec- 
trons are sorted by their increasing eigenenergies Sn^a- 
This ordering leads to a Liouvillian with a compact block 
structure. 

Introducing the abbreviation paa = {ci\p\o) for the el- 
ements of the density matrix in the energy eigenbasis 
I a) := |n, a) allows to write down the effective rate equa- 
tion for the populations [ ] in the form 



7ab,abPbb 



^ ^ ^ba,ba 
L b 



(A2) 



The transition rates from state b to state a are defined 
via 

8 

lab,ab= ^ la(3{£b-ea){a\Aa\b){a\A'^^\by , (A3) 

a, 13=1 

with transition matrix elements involving the system op- 
erators 

Al = di^L, A2 = d\^-^, As = di^R, A4 = d\ j^, 

is = d2,L, Ae = 4,L^ = d2,R^ As = dl j^. (A4) 

The coefficients 7^/3 (cj) in Eq. (A3) correspond to the 
even Fourier transformed bath correlation functions 



J —00 



(A5) 



which involve the bath operators 

k k 
k k 

k k 

= ^ h^R^k Ql^R^k^ ^8 = ^ ^2,i?,/e Q2,R,k' (A6) 
k k 

From this definition follows that the non- vanishing coef- 
ficients are given by 



712 (^ 
721 (^ 

734 (^ 
743 (^ 
756 (^ 
765 (^ 
778 {uJ 
787 (^ 



ri,L(-^)/i,L(-^), 

Fi,^ {uj) [1 - (uj)] , 
Ti^r{-u) fi^R{-u) , 
Ti^r{uj) [l-h,R{uj)], 

^2,L (-^) /2,L (-^) , 
r2,L {OJ) [1 - /2,L {OJ)] , 

r2,i?(-^)/2,i?(-^) , 

r2,R {^) [1 - /2,i? (a;)] . (A7) 
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